Universality classes in anisotropic non-equilibrium growth models 
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We study the effect of generic spatial anisotropics on the scaling behavior in the Kardar-Parisi- 
Zhang equation. In contrast to its "conserved" variants, anisotropic perturbations are found to be 
relevant in d > 2 dimensions, leading to rich phenomena that include novel universality classes and 
the possibility of first-order phase transitions and multicritical behavior. These results question 
the presumed scaling universality in the strong-coupling rough phase, and shed further light on the 
connection with generalized driven diffusive systems. 



PACS numbers: 05.40.+j, 64.60.Ak, 64.60.Ht. 

In non-equilibrium systems with conserved order pa- 
rameter, spatial anisotropics emerge as highly relevant 
perturbations, leading to drastic changes in the univer- 
sal scaling laws. For example, in driven diffusive systems 
(DDS) such as the driven Ising lattice gas jlj, an exter- 
nal drive generates a steady-state current that explicitly 
singles out one spatial direction. This leads to charac- 
teristic singularities in the structure factor already in the 
disordered phase, implies anisotropic scaling at the phase 
transition, and generates rich ordered structures in the 
low-temperature phase [0. Related anisotropic scaling 
behavior even ensues near the critical point of an Ising 
or Heisenberg system with conserved order parameter dy- 
namics that is driven out of equilibrium simply by impos- 
ing different Langevin noise strengths in different spatial 
directions, thus violating detailed balance 

Yet, one might expect that in models with a non- 
conserved order parameter, the picture could rather look 
similar to the situation in equilibrium systems, where 
anisotropies have much less dramatic effects. In order to 
further explore this issue, we study the Kardar-Parisi- 
Zhang (KPZ) model for kinetic roughening (5|, which has 
become another prototypical model for generic scale in- 
variance far from equilibrium HQ. Its anisotropic gen- 
eralization reads 
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with relaxation constants v\\/x > 0, but no restrictions 
on the signs of the non-linearities A||/j_ that describe 
curvature-driven growth. We denote the dimensionali- 
ties of the longitudinal and transverse sectors with dn and 
dx, respectively, with dii + dx = d. 77 is a stochastic driv- 
ing force with Gaussian correlations determined by the 
second moment (r?(x' ,t')n(x,t)) = 2D6(x — x')8(t — t'). 
After a simple length rescaling one may of course either 
choose i/n = vx or |An| = |Ax|- 

A previous study of this model restricted to two spa- 
tial dimensions (d\\ = dj_ = 1) found the anisotropy to 
be irrelevant for A|| and both positive fa]. Interest- 
ingly, if the non-linear terms come with opposite sign, the 



anisotropy is still irrelevant, but in addition also the non- 
linearity itself scales to zero. In a related model, where 
one of the non- linear couplings vanishes, the relevance of 
the remaining non-linearity was shown to depend on the 
values of dn and d± j|. It would therefore appear that 
anisotropies play only a minor role in kinetic roughening 
phenomena: Only if we render the signs of the non- linear 
couplings different, or set one of them to zero, do we seem 
to obtain any change in the scaling behavior. We shall, 
however, see that in any physical dimension larger than 
2, the above anisotropies destroy isotropic scaling at long 
wavelengths and in the long-time limit, and generate re- 
markably diverse behavior. This supports earlier conjec- 
tures that perhaps the strong-coupling scaling regime in 
the KPZ problem is not governed by universal scaling ex- 
ponents [fTo|], but rather subtly depends on microscopic 
details. These are of course captured only rudimcnta- 
rily in our simple anisotropies. For the original KPZ 
non-equilibrium growth or driven interface problem, this 
issue would appear to be mostly of academic interest; yet 
in light of a recent suggestion that the asymptotic scal- 
ing properties of interfaces pulled into unstable regions 
should be described by the (d + l)-dimensional rather 
than the d-dimensional KPZ equation [jy], our findings 
may well become directly accessible to real experiments. 

We have analyzed the non-linear Langevin equation 
([!]) by means of the dynamic renormalization group to 
one-loop order in the perturbation expansion with re- 
spect to the non-linear couplings A||/j_. The calculation 
is a straightforward generalization of the one-loop treat- 
ment for the isotropic KPZ equation or the equivalent 
noisy Burgers equation ^J^]. We map the Langevin 
equation to a dynamic functional, and proceed therefrom 
using standard field-theoretic tools The renormal- 

ization constants that track the ultraviolet singularities 
in d > 2 dimensions arc determined at a finite exter- 
nal wave vector q = ji (or cquivalently, at non-zero ex- 
ternal frequency) |L4|]. The scaling behavior of the the- 
ory with respect to this normalization scale is then en- 
coded in the associated RG equations, which are solved 
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by the method of characteristics (/i — > /iZ). The ensuing 
first-order differential flow equations define the running 
couplings as functions of the flow parameter / (or mo- 
mentum scale /i). We have to distinguish between the 
cases where both non-linearities have the same or oppo- 
site signs, A = Am = ±Aj_. The RG recursion relations 
Ida/dl — I3 a (g, 7) for the anisotropy ratio 7 = v\\/v±_ and 
the effective coupling g = AdfJ, d ~ 2 \ 2 D^f d± / 2 where 

A d = T(2 - ±)/[2 d+1 ir d /' 2 d{d + 2)] denotes a geometric 
factor, are determined by the RG beta functions 



Pg =g[d-2 + ( D - 3G,„ 



with the explicit one-loop flow functions 

Co = -g h (d|| + 2) + d x (d_L + 2) 7 2 ± 2d\\dxi] 
Cr/y = 3 [d(d|| ± d±7) - 4] , 
Go. = 91 [d(±d\\ + d ± j) - 4 7 ] . 



(2) 
(3) 

(4) 
(5) 
(6) 



The recursion relations for the special case A^ = arc 
readily obtained by setting 7 = in the expressions for 
these zeta functions. The zeros of the beta functions (j|) , 
(||) yield the RG fixed points, describing scale-invariant 
behavior. The universal infrared scaling exponents are 
then given by the corresponding fixed-point values of 
the zeta functions (Q)-(^|). Below the critical dimension 
d c = 2, the theory is ultraviolet-finite, and the emerg- 
ing RG fixed point is infrared-stable. For d > d c , the 
field theory remains renormalizable as a consequence of 
the emergence of an infrared- unstable fixed point and 
the underlying infinitesimal tilt invariance of the inter- 
face (Galilean invariance for the Burgers equation) [ p"5| |. 
Technically, renormalizability applies only in a system- 
atic e expansion about the lower critical dimension for 
the roughening transition (e = d— 2) . We shall neverthe- 
less consider the theory also at fixed d = dt\ + d±, with 
either sector dimension or A = d\\ — d± as parameters. 

The first question to be addressed is whether the 
anisotropy is a relevant perturbation at the isotropic 
KPZ fixed point. A previous study established that 
the isotropic fixed point is stable in 2 + 1 dimensions 
(d\\ = d± — 1) ||. However, our one-loop analysis re- 
veals this to be far from true in general. In fact, the RG 
flow equations, specifically Eq. (^) with (|^) and (ph, al- 
low for four different scenarios, depending on the values 
of d and A. Without loss of generality we can restrict 
ourselves to A > 0; for A < merely the roles of per- 
pendicular and parallel components are switched. In an 
analysis at fixed d and A we find the following regimes 
(depicted in Fig. Q): 

In regime A, d < -A/2 + ^(A/2) 2 + 8, both the 
isotropic fixed point 71 = 1 and the anisotropic fixed 
point ja = (4 — dd\\)/(dd± — 4) < are stable with an 
intermediate unstable uniaxial fixed point 70 = 0. Note 



(see remark after flow equations) that a negative fixed 
point in the anisotropy ratio does not mean that one of 
the relaxation constants v a becomes negative, but that 
we are at a fixed point where the non-linearities have op- 
posite sign. For the special case d\\ = d± one obtains 
7a = —1, and thus Wolf's earlier results || are recov- 
ered. In regime B, -A/2 + v/(A/2) 2 + 8 < d < V8, the 
anisotropic fixed point becomes positive, < ja < 1, 
and switches roles with the uniaxial fixed point 70. In 
regime C, y/8 < d < A/2 + ^(A/2) 2 + 8, the anisotropic 
fixed point becomes 7^4 > 1 and again stable, whereas the 
isotropic fixed point is now unstable. Finally, in regime 
D, d > A/2 + ^/(A/2) 2 + 8 , the anisotropic fixed point 
assumes again a negative value. Yet here the stability 
features of the fixed points are the converse of regime A. 
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FIG. 1. Stability diagrams for the fixed points of the 
anisotropy ratio 7 as functions of d and A; also indicated 
are (d\\,d±) at some particulary interesting dimensions. 

Note that the stability boundaries for the various fixed 
points of the anisotropy ratio 7, as determined from 
Eq. (^), merely require the existence of a finite fixed point 
< g* < 00, and are independent of its actual value. 
Hence one may speculate that a topologically similar 
if not identical stability diagram applies for the critical 
fixed point g c even beyond our one-loop approximation, 
and might perhaps even extend to the "strong-coupling" 
scaling behavior in the rough phase. 

The RG flows, projected onto the 7-g plane, are illus- 
trated in Fig. ||for scenario C. There are four critical fixed 
points, each one representing a different universality class 
(as discussed below). In the neighborhood of the uniax- 
ial and the anisotropic fixed points, the flows along the 
critical surface are towards each of them, since they each 
possess only one relevant scaling variable, the effective 
coupling constant g. The fixed point at 71 = 1 (describ- 
ing the isotropic kinetic roughening transition), located 
at the junction of two lines of critical points, is an exam- 
ple of a non- equilibrium bicritical point. The RG flows 
are quite reminiscent of a bicritical point in uniaxial mag- 
netic systems, with the high- and low-temperature phases 
of the magnetic system corresponding to the smooth and 
rough phases of the kinetic roughening problem. By anal- 
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ogy this suggests that there might actually exist a first- 
order phase transition along the line 71 = 1 for g > g c , 
separating two distinct non-equilibrium strong-coupling 
phases of different symmetry. Whereas the latter state- 
ment is somewhat speculative, one can systematically 
study the critical behavior of the new universality classes 
at the roughening transition. 

If 7 — > 0, i.e. V\_ — ► 00, the model reduces to the 
uniaxial KPZ equation with = 0. This can be seen 
by direct inspection of the RG flow equation, but it is 
also suggested by intuition since an infinitely large sur- 
face tension v±_ will always dominate over any finite non- 
linearity. With u = d\\h one realizes immediately that 
this uniaxial KPZ equation maps onto the generalized 
driven diffusion equation (GDDS) 

a t u = (^a2 + ^ai)u + ^ia||u 2 + c, (7) 

where u is a vector field with d\\ components and £ is 
a conserved noise with correlator (Ci(x, t)Q(x.' ,t')) = 
-2DS i:i V 2 d(x - x')<S(i - t'). Eq. (0) was originally in- 
troduced and analyzed to describe driven line liquids M . 
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FIG. 2. Sketch of the projected renormalization group 
flows for the anisotropic KPZ equation for scenario C. 

A crucial observation is that the nature of the fixed 
point of the effective coupling constant g changes dras- 
tically as a function of d and A. This can be seen by 
inspecting the RG beta function for g in Eq. (||). Note 
that the fixed-point value of the effective coupling con- 
stant diverges and changes sign at Cd = (3 — d±/2)^. 
In an e expansion with respect to the critical dimension 
d c = 2 one then finds two completely different regimes: 
For Co < (3 — dx/2)Ci/N) i- e -> to leading order in e for 
A > 0.47 — 1.89e, there exists a critical fixed point of 
order e = d — 2 above two dimensions marking a non- 
equilibrium phase transition from a smooth to a rough 
phase. In the converse case, the non-linearity is pertur- 
batively irrelevant above d c = 2, and consequently the 
roughening transition disappears. In addition, below two 
dimensions one finds a stable fixed of order e' = 2 — d 
describing anomalous scaling behavior. In other words, 
there exists a critical value A c « 0.47 where d c = 2 
changes its role from a lower to an upper critical dimen- 



sion, and the elusive strong-coupling scenario is trans- 
formed into a weak-coupling one. 

These observations have several implications. First, 
owing to the GDDS connection, it is possible to access 
the one-dimensional KPZ scaling (du = 1, d± = 0) per- 
turbatively by expanding around the upper critical di- 
mension d c — 2 of the standard non-critical DDS model 
{d\\ = 1, d± — I). Second, a naive extrapolation of 
A = 0.47- 1.89e toe = 1 (d = 3) would suggest that DDS 
with cZ|| = I and d± = 2 might display a non-equilibrium 
phase transition. However, this is most likely an artifact 
of the e expansion, and higher-order terms in the pertur- 
bation series will probably result in a boundary line for 
the critical fixed point that asymptotically approaches 
d\\ = I as d± -> oo such that purely mean-field behavior 
ensues for d > 2. Yet careful numerical simulations of 
the GDDS model that map out the (d, A)-plane and de- 
termine the exact location of the boundary line, beyond 
the regime close to d — 2 accessible in a perturbative RG 
approach, would certainly be very desirable. 

There are now several ways to calculate the expo- 
nents in an e expansion with respect to d c = 2. Keep- 
ing, e.g., A = dn — d± fixed, we find z\\ = 2 + 
2 (A - 2) e/(A 2 + 8A - 4) and z± = 2, which describes, 
depending on the value of A, the critical fixed point of 
the roughening transition or the stable fixed point below 
d c . For example, extrapolating to (dy, d±) = (2, 1) yields 
Z|| 1.6. The actual value for z found here is actually 
not crucial — the one-loop analysis should not be ex- 
pected to be very accurate — but the very fact that it 
differs from z c = 2, the exact result for the roughening 
transition of the isotropic KPZ equation. Furthermore, 
the RG analysis at the uniaxial fixed point shows no sign 
of an upper critical dimension, while it is known that 
d — 4 is the upper critical dimension for the roughening 
transition of the isotropic KPZ model H , E6[ . 

Besides the uniaxial fixed points 70 = and 7^ = 00 
there exists also an anisotropic fixed point with a fi- 
nite anisotropy ratio ja = (4 — ddu )/(ddx — 4). We 
find that the non-equilibrium dynamics described by the 
anisotropic and isotropic fixed point are markedly differ- 
ent. At first sight this may seem quite surprising since 
at both fixed points the parallel and perpendicular sur- 
face tension term v\\ t ± show identical scaling behavior 
under a RG transformation; the only difference resides 
in their amplitudes. The reason for this unusual behav- 
ior is that 7 = 1 is a non- generic high-symmetry point in 
parameter space, allowing the isotropic KPZ problem to 
be mapped onto the statistical mechanics of a directed 
polymer in a random medium. The high symmetry is 
also reflected in the mean-field value for the dynamic ex- 
ponent z c = 2 at the non-trivial fixed point describing 
the roughening transition. In contrast, our RG analysis 
at the anisotropic fixed point yields a dynamic critical 
exponent different from 2, to leading order in a 2 + e ex- 
pansion at fixed A given by z = 2 + e (A 2 -4)/(10A 2 -8). 
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Similar to the GDDS equation, there exists a critical line 
A c (d) such that for |A| < A c (d) there is no roughening 
transition above d = 2. To leading order in e we ob- 
tain A c {d) = ±(2/V5- 13e/5\/5). Since \A c (d)\ < 1 
for d > 2, a scenario with no roughening transition ap- 
pears unlikely for any integer dimensions d|| and d±. 
For A > A c (d) there is, as for the uniaxial fixed point, 
one set of dimensions which is particularly interesting, 
namely (dii,dx) = (2, 1). Upon extrapolating the result 
obtained from our e expansion we arrive at an unphysi- 
cally low critical value z c = 1/2 ± 0(e 2 ). Unfortunately, 
this leading-order result is not accurate enough to wager 
a reliable quantitative prediction for the actual value of 
the dynamic critical exponent. But, as before, the im- 
portant conclusion to be drawn from the RG analysis is 
again that the dynamic critical exponent differs from the 
isotropic value z c — 2. In order to gain reliable quantita- 
tive estimates for the exponents of this novel universality 
class, careful numerical simulations would be invaluable. 

From the above analysis of the various universality 
classes and our understanding of the isotropic KPZ equa- 
tion we derive the following picture. With the exception 
of the particular case d\\ = dj_, the anisotropic fixed point 
always displays strong-coupling behavior, i.e. we find a 
roughening transition above the lower critical dimension 
d c = 2 similar to the isotropic fixed point. However, the 
universality class of the roughening transition is markedly 
different from the isotropic KPZ model. In contrast to 
the anisotropic and the isotropic fixed points, the uniax- 
ial fixed point 70 = (GDDS) shows extended regions in 
the (d, A) plane with strong- and weak-coupling behav- 
ior, respectively. In particular, for d|| < d± the roughen- 
ing transition is entirely absent and we have anomalous 
scaling only below d c = 2. 

Some specific examples serve to illustrate the impli- 
cations of the above results. For (d\\,d±) = (3,1) the 
isotropic fixed point 71 = 1 is unstable. If 7 < 1, the RG 
flow tends towards the uniaxial fixed point. The corre- 
sponding GDDS universality class lacks the special sym- 
metry of the isotropic KPZ equation implying that the 
dynamic critical exponent differs from z c = 2. If 7 > 1, 
the RG trajectories flow towards 7 = 00. This fixed point 
lies again in the GDDS universality class, but with inter- 
changed roles of (iy and d±. Because of the symmetry of 
the model, the asymptotic behavior is equivalent to that 
at the fixed point 70 = in (du,d±) = (1,3) space di- 
mensions, where no roughening transition occurs above 
d c = 2. This implies not only that the isotropic KPZ 
equation is unstable in d = 4 with respect to even the 
slightest amount of anisotropy, but more dramatically 
also that there exists no phase transition and hence no 
strong-coupling behavior at all, provided the relaxation 
amplitudes satisfy V\\ < for d\\ =1. Similar critical 
end points are found throughout domain D of Fig.[l]. 

Finally, we note that we have performed a similar one- 
loop RG analysis for the "conserved" KPZ variants, with 



essentially an additional Laplacian on the right-hand-side 
of Eq. (Q), both with conserved JlTj and non-conserved 
white noise Neither of these models displays a 

roughening transition. Rather, there is a single universal 
scaling regime with non-trivial exponents below the up- 
per critical dimensions d c = 4 and 2, respectively. It then 
turns out that anisotropic perturbations do not matter 
crucially: The standard isotropic fixed point remains sta- 
ble in both cases [[l9| . It would thus appear that the novel 
features reported above for the anisotropic KPZ equa- 
tion are crucially connected with the very existence of a 
roughening transition, and thus with the scaling proper- 
ties of a strong-coupling rough phase. 
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